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ABSTRACT 

Supernova remnants may exhibit both thermal and nonthermal X-ray emission. Such rem- 
nants can be distinguished by the weakness of their X-ray lines, because of the presence of a 
strong nonthermal X-ray continuum. RCW 86 is a remnant with weak lines, resulting in low and 
peculiar abundances when thermal models alone are used to interpret its X-ray spectrum. This 
indicates the presence of a strong nonthermal synchrotron continuum. 

We analyze ASCA X-ray spectra of RCW 86 with the help of both nonequilibrium ionization 
thermal models and nonthermal synchrotron models. A two-temperature thermal model and a 
simple nonthermal model with an exponential cutoff (plus interstellar absorption) give reasonable 
results. We obtain blast wave velocity of 800 km s _1 , the shock ionization age of 1 — 3 x 10 11 
cm~ 3 s, and the break in nonthermal spectra at 2 — 4 x 10 16 Hz. The strength of nonthermal 
continuum correlates well with the radio brightness in the bright SW section of the remnant. 
This is convincing evidence for X-ray synchrotron emission in RCW 86. 



Subject headings: ISM: individual (RCW f 

RCW 86 AND ITS X-RAY SPEC- 
TRUM 



RCW 86 (G315.4-2.3) is a shell-like supernova 
remnant (SNR) in the southern sky, 42' in angu- 
lar size, and 35 pc in linear size at a kinematic 
distance of 3 kpc found from optical observations 
by Rosado et al. (1996). (The distance estimate is 
uncertain, and values as low as 1 kpc cannot be ex- 
cluded according to Smith 1997). Optical emission 
from radiative shocks is present at various loca- 
tions along its shell, with by far the brightest com- 



supernova remnants - X-rays: ISM 

plex of optical filaments in the SW. This bright 
complex, to which the optical designation "RCW 
86" actually refers, has been generally interpreted 
as the result of the interaction of the blast wave 
with a dense (~ 10 cm -3 ) interstellar cloud (e. g., 
Rosado et al. 1996; Smith 1997). The slow (~ 100 
km s _1 ) radiative shocks are thought to have been 
driven into the cloud by the blast wave impact. 
Optical, Balmcr-dominated emission from faster 
nonradiative shocks was also detected (Long & 
Blair 1990), and Smith (1997) found that Balmer- 
dominated optical filaments almost completely en- 
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Table 1: 1997 ASCA Observations of RCW 86 



Name 


RA (2000) 


DEC (2000) 


Sequence 


Time 


Duration 


Effective SIS 




(hh mm ss) 


(° ' ") 


Number 


(yy.ddd) 


(ksec) 


Exposure (ksec) 


RCW 86 (FRONT) 


14 40 34.8 


-62 41 17 


55045000 


97.048 


136 


39 


RCW 86 (BACK) 


14 40 54.9 


-62 38 48 


55046000 


97.050 


121 


42 


RCW 86 (BK SI) 


14 40 18.4 


-62 42 24 


55046010 


97.226 


134 


44 



circle RCW 86. This is presumably a blast wave 
propagating into a low-density (<~ 0.2 cm~ 3 , Long 
& Blair 1990) interstellar medium (ISM). Its ve- 
locity is in the range of 400-900 km s -1 , found 
from the width of the broad Ha component (Long 
& Blair 1990; Ghavamian 1999; Ghavamian et al. 
2000). The remnant's age is ~ 10 4 yr, derived 
from the remnant's angular size of 42', its distance 
of 3 kpc and the known shock speed (Rosado et 
al. 1996). (However, note that RCW 86 may be 
as young as 2000 yr if it is located at distance 
of 1 kpc.) As expected, thermal X-ray emission 
produced behind fast shocks (Pisarski, Hclfand, & 
Kahn 1984), nonthermal radio synchrotron emis- 
sion from relativistic electrons (Whiteoak & Green 
1996: Molonglo catalog), and strong IR emission 
from collisionally heated dust (Dwek et al. 1987; 
Arendt 1989; Greidanus & Strom 1990; Saken, 
Fesen, & Shull 1992) are all present. RCW 86 
was thought to be a remnant of SN 185 (Clark 
& Stephenson 1977), but even the existence of 
this SN has been questioned (Chin & Huang 1994; 
Schaefer 1995). 

RCW 86 is clearly an interesting SNR, seen 
across the wavelength band from radio to X-rays, 
with both slow and fast shocks present. But the 
most interesting is its peculiar X-ray spectrum 
as revealed by the ASCA Performance Verifica- 
tion (PV) observations (Vink, Kaastra, & Blocker 
1997). The X-ray lines are generally much weaker 
than expected from a normal (solar) abundance 
plasma, with the exception of the Fe Ka line 
at 6.4 keV. Vink et al. (1997) fitted ASCA 
data with two-temperature nonequlibrium ioniza- 
tion (NEI) models and found strongly subsolar 
abundances (except for Fe in their hot tempera- 
ture component) and extremely low (< 10 10 cm~ 3 
s) ionization timescales n e t. These results are 
clearly unphysical, and Vink at al. speculated 
that perhaps they are caused by strong devia- 



tions from the Maxwellian electron distribution for 
low-temperature plasmas. A more recent study 
of RCW 86 by Bocchino ct al. (2000), based on 
BcppoSAX satellite data, confirmed the weakness 
of emission lines (other than the Fe Ka line) , al- 
though the authors attributed this effect to ele- 
mental depletion onto dust. 

The goal of this paper is to explain the puzzling 
X-ray spectrum of RCW 86, based on new, high 
signal-to- noise (S/N) ASCA observations. These 
new observations are presented and briefly dis- 
cussed in § 2, while in § 3 we analyze these obser- 
vations in detail and show that they can be succes- 
fully explained as a combination of thermal emis- 
sion and nonthermal synchrotron emission. In § 4, 
we summarize our results and discuss prospects 
for further progress. 

2. 1997 ASCA OBSERVATIONS OF RCW 
86 AND THEIR INTERPRETATION 

Three long observations of RCW 86 were per- 
formed by ASCA in 1997 (Tabic 1). The data 
were processed under standard Revision 2 Data 
Processing (REV2) and screened according to the 
following criteria: elevation angle greater than 10 
degrees, geomagnetic field structure with cut-off 
rigidity greater than 6, observations outside South 
Atlantic Anomaly, stable pointing, and stable in- 
ternal background (from monitored Gas Imaging 
Spectrometers count rates). The effective expo- 
sure times after screening are typically about 40 
ks as shown in Table 1. Note that the effective ex- 
posure time of the PV observations used by Vink 
ct al. is only 14 ks. The long exposure times of 
the new observations resulted in good S/N spec- 
tra, leading to a significant improvement over the 
PV observations, particularly at intermediate and 
high photon energies. This allowed us to separate 
thermal and nonthermal contributions to the X- 
ray spectrum. 
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Fig. 1.— ROSAT PSPC image of RCW 86. Fields 
of view of two ASCA observations from 1997 (SIS 
- squares, GIS - circles) cover the brightest, SW 
"knee" region of the remnant. 
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Fig. 2.— GIS image of the SW knee of RCW 86, 
overlayed with 0.843 GHz radio contours (based 
on MOST radio data - Whiteoak & Green 1996). 
The X-ray and radio emission peaks nearly over- 
lap, but radio emission is much stronger in the 
horizontal arm to the west of the knee bend than 
in the vertical arm to the north. 



Fields of view of solid-state imaging spectrom- 
eters (SIS) and gas-imaging spectrometers (GIS) 
for the first two observations listed in Table 1 (se- 
quences 55045000 and 55046010) are plotted on 
the ROSAT PSPC image of RCW 86 (Figure 1). 
The third pointing in the SW (sequence 55046000) 
is not shown, because its SIS field of view is cen- 
tered outside the remnant. However, this pointing 
contains useful GIS data on the SW part of the 
remnant. These three long pointings provide high 
S/N images and spectra of the brightest part of 
RCW 86. GIS and SIS images are shown sep- 
arately in Figures 2 and 3, respectively, while a 
GIS hardness map defined as the soft (< 3 keV) 
to hard (> 3 keV) count rate ratio is shown in 
Figure 4. The GIS hardness map is affected by 
the difference in the SIS and GIS point spread 
functions. Significant spectral variations are read- 
ily apparent in this map, in accord with findings 
of Vink et al. (1997) which were based on PV 
observations with a much lower S/N ratio. The 
bright "knee" in the SW corner of the remnant 
exhibits striking variations in the hardness ratio, 
with the soft vertical arm and the progressively 
harder horizontal arm as one moves away (to the 
west) from the "knee" bend. The 0.843 GHz radio 
contours shown in Figures 2 and 3, from the Mo- 
longlo Observatory Synthesis Telescope (MOST) 
catalog (Whiteoak & Green 1996), demonstrate 
that the X-ray hardness is positively correlated 
with the radio/X-ray brightness ratio. 

SIS spectra (Figure 5) of three rectangular re- 
gions shown in Figure 4 demonstrate spatial vari- 
ations quantitatively. Spectra of the vertical and 
horizontal arms, which we label the Soft and 
the Hard Regions, respectively, differ dramati- 
cally. The Soft Region is brighter at low energies, 
while the situation reverses at high energies. The 
Medium-Hard Region at the knee bend is interme- 
diate in its hardness. Lines from O, Ne, and Mg 
may be discerned in these spectra at low energies, 
although they are apparently weak. Their pres- 
ence indicates that X-ray spectrum of RCW 86 is 
at least partially of the thermal origin. 

We first test the hypothesis that X-ray spec- 
tra in Figure 5 are produced entirely by thermal 
plasma, using a set of noncqulibrium ionization 
(NEI) models developed by us and publicly avail- 
able in the XSPEC vll software package. (For in- 
troduction to XSPEC, sec Arnaud 1996 or visit the 
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Fig. 3.— SIS image of the SW knee of RCW 86, 
overlayed with 0.843 GHz radio contours. The 
bright, nearly horizontal ridge of radio emission 
is located interior of the brightest X-ray emission. 
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Fig. 4.— The GIS hardness map, soft (< 3 keV) 
to hard (> 3 keV) count rate ratio, with soft re- 
gions in black and hard regions in white, overlayed 
with GIS intensity contours. Three SIS extraction 
regions are shown (clockwise, from the top: Soft 
Region, Medium-Hard Region, and Hard Region). 
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Fig. 5. — SIS spectra for the three extraction re- 
gions (Medium-Hard - top, Soft - middle, Hard - 
bottom, as seen below 1 keV) . Weak emission lines 
are clearly seen in all spectra. 



web site at http://heasarc.gsfc.nasa.gov.) Ther- 
mal NEI models are based on an updated Hamil- 
ton & Sarazin (1984) spectral code, which includes 
recent atomic calculations for Fe L-shell lines by 
Liedahl, Osterheld, & Goldstein (1995). These 
models are discussed in more detail by Borkowski, 
Lyerly, & Reynolds (in preparation). We used Se- 
dov models with the variable metal abundances 
(but with their relative abundances fixed in solar 
ratios) to fit these spectra, with the best fit mod- 
els shown in Figure 6. Shock temperatures, post- 
shock electron temperatures, and SNR ionization 
timescales (equal to the product of postshock elec- 
tron density and the remnant's age) are: 1.4, 1.2, 
and 2 keV; 1.3, 1.2, and 0.7 keV; and 5.0 x 10 10 
cur 3 s, 5.4 x 10 10 cm" 3 s, and 8.9 x 10 10 cm" 3 
s, for the Hard, Medium-Hard, and Soft Regions, 
respectively. In all three regions we arrived at the 
absorbing column number density Nh of 3 x 10 21 
cm~ 2 . Metal abundances in these models are low 
(0.2, 0.4, and 0.55 solar for the Hard, Medium- 
Hard, and Soft Regions, respectively), and spa- 
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Fig. 6. — Fits to SIS spectra with the Sedov 
models, shown as solid lines (Hard Region - top, 
Medium-Hard - middle, Soft - bottom). The spec- 
trum of the Hard Region was shifted upward by 
factor of 10 for the clarity of the presentation. 
Metal abundances in Sedov models are low (0.2- 
0.55 solar) and increasing from top to bottom. 



tially varying. Spatial abundance variations at 
the periphery of RCW 86 are not expected, and 
something is obviously wrong with these mod- 
els. Fits with thermal models with variable abun- 
dances presented by Vink et al. (1997) also lead 
to strongly subsolar abundances of heavy elements 
and extremely low ionization timescales. As these 
authors pointed out, this is clearly unphysical, and 
they also concluded that standard thermal mod- 
els do not provide a satisfactory description of the 
X-ray spectrum of RCW 86. 

The weakness of spectral lines suggests the 
presence of nonthermal X-ray continuum in RCW 
86, as has been recently inferred for SN 1006 
(Koyama et al. 1995), Cas A (Allen ct al. 1997), 
G347.5-0.2 (Slane et al. 1999), and possibly even 
in RCW 86 itself based on its high-energy spec- 
trum obtained by the Rossi X-ray Timing Explorer 
(Allen, Gotthelf, & Petre 1999). The spatial corre- 



lation between the radio and X-ray emission noted 
above suggests that this continuum may be pro- 
duced by synchrotron X-ray emission from rela- 
tivistic electrons. In order to test this hypothe- 
sis, we modeled SIS spectra as a sum of nonther- 
mal and thermal emission (plus interstellar ab- 
sorption). We used a simple nonthermal model 
with a power-law electron energy distribution and 
an exponential cutoff at high electron energies (the 
srcut model in XSPEC vll, see Reynolds 1998 and 
Reynolds & Keohane 1999 for more details) and 
a thermal Sedov model. The srcut model is pa- 
rameterized by a 1 GHz flux density (in Jy) and 
spectral index, and a break frequency which is the 
peak frequency emitted by electrons with the e- 
folding energy of the exponential cutoff (at which 
the spectrum has dropped by a factor of about 
3.7). Abundances were assumed to be 2/3 solar as 
appropriate for the ISM in the solar neighbour- 
hood (e. g., Snow & Witt 1996; Mathis 1996). 
Nonthermal models were normalized by the 1 GHz 
radio flux obtained by integrating 0.843 GHz ra- 
dio emission over the SIS extraction regions and 
extrapolating to 1 GHz using the observed radio 
spectral index of 0.6 (Green 1998). However, for 
the fitting, the radio-to- X-ray spectral index was 
allowed to vary. Resulting fits are shown in Fig- 
ure 7. The quality of the fits is better than for the 
pure thermal models, although they are still not 
statistically acceptable. But unlike for pure ther- 
mal fits, parameter values of these mixed thermal 
and nonthermal models are physically reasonable. 
The Sedov model parameters are very similar in 
all three regions: shock temperatures in the range 
from 1.0 to 1.4 keV, low (< 0.2 keV) postshock 
electron temperatures, and ionization timescales 
in the range from 1.4 x 10 11 cm~ 3 s to 1.8 x 10 11 
cm -3 s. For the nonthermal srcut models we ob- 
tain break frequencies of 7.5 x 10 16 , 8.2 x 10 16 , 
2.5 x 10 17 Hz for the Hard, Medium-Hard, and 
Soft Regions, respectively, and radio spectral fre- 
quencies in the range from 0.62 to 0.67, in sur- 
prisingly good agreement with the (uncertain) ob- 
served RCW 86 spectral index of 0.6. The lower 
electron temperatures are a natural result of at- 
tributing some of the continuum to nonthermal 
emission, but we stress that these models are not 
well constrained and this may be an artifact. 

But the actual situation is more complex than 
described above, because a weak Fe Ka line at 6.4 
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Fig. 7. — Fits to SIS spectra with the sum of Sedov 
models and srcut nonthermal synchrotron mod- 
els. Nonthermal contribution, shown separately 
for each region, is most important in the Hard re- 
gion (top), and least important in the Soft region 
(bottom). 



keV seen in the SIS spectrum of the Medium-Hard 
Region is not present in the models. This line, al- 
ready noted by Vink et al., is most likely produced 
by a high-temperature thermal plasma. In order 
to investigate this line in more detail, we combined 
available GIS2 data to construct a high S/N spec- 
trum (Figure 8), which clearly shows Fe Ka line. 
At high energies above 5.5 keV, these data are 
fit well by a bremsstrahlung continuum with elec- 
tron temperature kT e = 5.3 keV (with 3.7-8.8 keV 
90% confidence interval) plus a narrow Gaussian 
line at 6.42 keV (with 6.35-6.45 keV 90% confi- 
dence interval). Such a low energy of the Fe Ka 
line suggests a very low ionization timescale of the 
hot plasma. This is confirmed by fitting a sin- 
gle ionization timescale, single temperature ther- 
mal NEI model (vnei model in XSPEC, with vari- 
able Fe abundance and solar abundances for other 
heavy elements) to our GIS data above 5.5 keV. 
We obtain plasma temperature of 5.3 keV (with 



3.6-8.6 keV 90% confidence interval), a very low 
ionization timescale of 5.2 x 10 8 cm~ 3 s (less than 
6.7 x 10 9 cm~ 3 s with 90% confidence), and Fe 
abundance of 3.5 solar (with 1.2-3000 90% confi- 
dence interval). Model parameters and confidence 
intervals depend quite sensitively on the chosen 
energy range, resulting in progressively lower elec- 
tron temperatures as lower energy photons arc in- 
cluded in the fits. Atomic data for Fe Ka line 
emission produced by electron excitation of low- 
ionization stages of Fe are poor, further contribut- 
ing to uncertainties in derived parameter values, 
the Fe abundance in particular. It is also possible 
that Fe Ka is produced in dust grains (Borkowski 
& Szymkowiak 1997), which would make our es- 
timates of the ionization timescale unreliable. In 
this situation, the true confidence intervals for all 
model parameters are certainly much larger that 
the confidence intervals just quoted. In particular, 
approximately solar abundance of Fe is consistent 
with the current data. 

The best-fit thermal bremsstrahlung and nei 
models produce similar spectra across the GIS 
spectral range (Figure 8), without any prominent 
emission lines except for the fluorescent Fe Ka 
line. In the nei model this is caused by an ex- 
tremely low ionization state of the gas, where all 
heavy elements have not yet been ionized to their 
H- and He-like stages (in the case of Fe even M- 
shell electrons must still be present because of the 
apparent absence of the Fe L-shell emission). The 
presence or absence of lines depends very sensi- 
tively on the ionization timescale of the gas. If we 
assume that the ionization timescale is equal to 
6.7 x 10 9 cm~ 3 s (the upper 90% confidence limit 
in our fits), and fit the nei model to the GIS2 data 
above 5.5 keV, the lines are already very promi- 
nent (Figure 8). This extreme sensitivity to the 
ionization timescale means that the contribution 
of the high-energy thermal component to the X- 
ray spectrum at low photon energies is not well 
constrained by the data, although in all cases this 
contribution is not negligible and must be included 
in spectral modeling. 

The presence of a poorly-constrained high- 
temperature thermal component in RCW 86 sig- 
nificantly complicates analysis of its X-ray spec- 
trum, and our previous arguments about the ne- 
cessity of the nonthermal component need to be re- 
examined. Is it possible that a mixture of gas with 
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Fig. 8. — The Fe Kq line at 6.4 keV is clearly seen 
in the high S/N GIS2 spectrum of the SW section 
of RCW 86. Three models providing acceptable 
fits to energies above 5.5 keV are shown: a best- 
fit single ionization timescale model nei with ion- 
ization timescale r of 5 x 10 8 cm~ 3 s (solid line), 
an nei model with r = 6.7 x 10 9 cm~ 3 s (dashed 
line), and a simple bremsstrahlung model with a 
Gaussian line at 6.4 keV (dotted line). 



standard ISM abundances but with widely differ- 
ent temperatures results in the observed weakness 
of the lines? The answer seems negative, according 
to Vink et al. (1997), and according to Bocchino 
ct al. (2000) who fitted BeppoSAX spectra of the 
SW region of RCW 86 with a two-component NEI 
model. Bocchino et al. (2000) derived the fol- 
lowing parameters for the high-temperature com- 
ponent: kT e = 5.7 keV, ionization timescale of 
1.9 x 10 8 cm~ 3 s, and the Fe abundance of 3.0 with 
respect to solar, in excellent agreement with our 
results. But the derived abundances in the low- 
temperature (1.0 keV) NEI component are again 
strongly (0.25) subsolar. We also attempted to 
fit two-component (Sedov and nei) thermal mod- 
els to our ASCA SIS spectra, but the the derived 
abundances in Sedov models are again strongly 



subsolar and spatially varying. Although this re- 
sult does not rule out the possibility of the pure 
thermal explanation, it demonstrates that such 
an alternative is unlikely, in accord with Vink 
ct al. (1997). Another alternative explanation 
for the apparent subsolar abundances, elemen- 
tal depletion onto dust (Bocchino et al. 2000), 
is not supported by the abundance patterns de- 
rived from ASCA and BcppoSax observations. In 
view of these findings, we now present our ef- 
forts to construct a physically reasonable model 
for the X-ray spectrum of RCW 86, in which we 
allow for a simultaneous presence of nonthermal 
synchrotron emission and thermal emission from 
multi-temperature plasma. 

3. X-RAY SPECTRA OF RCW 86 AS 
A MIXTURE OF THERMAL AND 
NONTHERMAL EMISSION 

We begin with the analysis and modeling of the 
brightest extraction region, the Medium-Hard Re- 
gion at the knee bend (see Figure 4). In addition 
to the SIS spectra described in the previous sec- 
tion, we also use the GIS spectrum extracted from 
the same spatial region in order to improve S/N 
ratio at higher photon energies. Below 1.2 keV, 
the PV SIS spectra are used instead of the newer 
SIS and GIS data because of their superior spec- 
tral resolution and good S/N ratio. These three 
data sets are shown in Figure 9. We always in- 
clude these data sets together in our analysis, by 
performing joint fits to the PV SIS data, and the 
1997 SIS and GIS data. In order to allow for mis- 
match between normalizations of individual data 
sets, we allowed them to vary when performing 
joint fits to the data. 

The thermal plasma emission is modeled as 
a sum of low-temperature and high-temperature 
components. For the low-temperature component, 
we use either a simple plane-parallel shock model 
with constant electron temperature kT e and ion- 
ization age r (the pshock model in XSPEC) or a 
much more complex spherically-symmetric shock 
with unequal ion and electron temperatures (the 
Sedov model) which was already described pre- 
viously. The high-temperature component is de- 
scribed by a simplest NEI model, a single ioniza- 
tion timescale model nei. The nonthermal syn- 
chrotron emission is modeled by the sreut model, 
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using the observed radio flux and its spectral in- 
dex, and allowing the break frequency to vary. (If 
a power law is used instead of the scrut model, the 
separation between the thermal and nonthermal 
components becomes more difficult, and the fitted 
power-law index and normalization of the power- 
law component are also difficult to interpret.) 

Because of the presence of synchrotron emis- 
sion, the chemical abundances cannot be reliably 
determined from the observed spectra, so that 
we assumed 2/3 solar abundances in the low- 
temperature thermal component. Furthermore, 
the S /N ratio of the current data sets at high pho- 
ton energies is not adequate for a unique separa- 
tion between the high-temperature thermal com- 
ponent and the nonthermal synchrotron emission. 
In this situation we assumed that all emission 
above 5.5 keV is thermal, and used parameters 
of the high-tcmpcrature thermal component ob- 
tained by fitting spatially-integrated GIS spectra 
above 5.5 keV (see § 2). The normalization of this 
component was obtained by fitting this thermal 
model alone to the spectral data above 5.5 keV. 
The modeling results for the two models described 
above are tabulated in Table 2. 

We generally obtain similar results for our 
two types of models, pshock+nei+sreut and Se- 
dov+nei+sreut, although fits are always better for 
the former. The best-fit pshock+nei+sreut model 
for the Medium-Hard Region is shown in Figure 9. 
The low- (high-) temperature thermal component 
dominates at low (high) energies, and the non- 
thermal synchrotron emission is most important 
in the intermediate energy range. However, the 
nonthermal component contributes at all energies, 
accounting for 22% of the total detected flux in the 
0.5 - 10 keV range (this fraction is shown in the 
last row of Table 2). The top panel in Figure 9, 
with the high-quality PV SIS data, shows most 
clearly how strong nonthermal emission leads to 
low equivalent widths for the prominent Mg and 
Si lines. 

Modeling of the Soft and Hard Regions is done 
in exactly the same way as for the Medium-Hard 
Region, except that we do not use PV SIS data 
because of their lower S/N. (For the Hard Region 
and the Sedov+nei+scrut model, the shock tem- 
perature and the postshock electron temperature 
in the Sedov model could not be determined inde- 
pendently, so that we fixed the postshock electron 
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Fig. 9. — ASCA spectrum of the Medium-Hard 
Region (PV SIS data - top, 1997 SIS data - mid- 
dle, GIS data - bottom). A three-component 
[pshock + nei + scrut) model is also plotted (total 
- thick solid line, thermal components - thin solid 
lines, nonthermal synchrotron component - dotted 
line. 



temperature to a low value of 0.05 keV, similar to 
the postshock electron temperatures found in the 
Medium-Hard and Soft Regions.) Our SIS and 
GIS data for the Soft Region, together with the 
best-fit pshock+nei+sreut model fom Table 2, are 
shown in Figure 10. As expected, the low- (high-) 
temperature thermal component dominates again 
at low (high) energies, but the nonthermal com- 
ponent contributes only 1.8% to the total flux in 
this model. 

The situation is very much different in the Hard 
Region (Figure 11), where nonthermal emission is 
very prominent, contributing 33% of the total flux. 
Again, emission at high energies is dominated by 
the high-temperature thermal component with an 
Fe Ka line clearly seen at 6.4 keV, and faint lines 
with low equivalent widths at low energies are pro- 
duced by the low-temperature component. But 
the emission measure of the latter is at least twice 
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Table 2: Spectral Models 



Medium-Hard Region Soft Region Hard Region 

Model Parameters A a B b A B A B 

N H (cm" 2 ) 1.9 x 10 21 2.5 x 10 21 2.8 x 10 21 3.3 x 10 21 3.2 x 10 21 3.3 x 10 21 

kTi (keV) 0.75 0.90 c (0.04 d ) 0.77 1.10 (0.07) 0.80 1.4 (0.05 e ) 

r ( (cm- 3 s) 1.8 x 10 11 3.1 x 10 11 1.1 x 10 11 2.2 x 10 11 5.7 x 10 10 1.5 x 10 11 

Abundance (fixed) 2/3 2/3 2/3 2/3 2/3 2/3 

EMi/(Aird 2 ) (cm" 5 ) 9.8 x 10 11 1.3 x 10 12 1.1 x 10 12 1.4 x 10 12 4.2 x 10 11 5.3 x 10 11 

kT h (keV) (fixed) 5.3 5.3 5.3 

T h (cm- 3 s) (fixed) 5.2 x 10 8 5.2 x 10 8 5.2 x 10 8 
Abundance (fixed) 111 

Fc (fixed) 3.4 3.4 3.4 

EM h /{Aird 2 ) (cm- 5 ) 2.6 x 10 11 1.2 x 10 11 2.4 x 10 11 

v c (Hz) 3.2 x 10 16 3.1 x 10 16 1.9 x 10 16 1.5 x 10 16 3.7 x 10 16 3.7 x 10 16 

a (fixed) 0.6 0.6 0.6 0.6 0.6 0.6 

Fighz (Jy) 4.48 4.48 0.57 0.57 3.19 3.19 

X f (0.5-10 keV) 22% 20% 1.8% 1.2% 33% 33% 



a pshock + nei + srcut 
b Sedov + nei + srcut 
c Postshock gas temperature 
d Postshock electron temperature 
e Fixed 

f Nonthermal/ (nonthermal+thermal) flux ratio 



as low as in the Medium-Hard and Soft Regions 
(see Table 2). At the same time, the nonthermal 
component appears to be at the same intensity 
level as in the Medium-Hard Region, leading to 
the increased fraction of the nonthermal emission 
in the Hard Region. 

The nonthermal synchrotron emission is strongly 
correlated with the radio emission. In the 
pshock+nei+sreut model, after correction for ab- 
sorption, we obtain the total 0.5-10 keV X-ray 
synchrotron fluxes of 9.5 x 10 -12 ergs cm -2 s _1 , 
6.1 x 10- 13 ergs cm -2 s _1 , and 8.1 x 10~ 12 ergs 
cm- 2 s _1 for the Medium-Hard, Soft, and the 
Hard Regions, respectively. The corresponding 
1 GHz radio fluxes are 4.5 Jy, 0.57 Jy, and 3.2 Jy. 
The low synchrotron X-ray flux is accompanied 
by the low radio flux in the Soft Region, while 
an order of magnitude increase in the nonthermal 
X-rays in the Medium-Hard and Hard Regions is 
matched by a similar increase in the radio flux. 
Note that such drastic changes are not seen in the 
thermal X-ray flux in the 0.5-10 keV range, which 



varies by only 20% between the regions. Such 
strong positive correlation between nonthermal X- 
ray and radio emission is seen in the prototypical 
synchrotron-dominated remnant of SN1006. The 
presence of this correlation in RCW 86 strongly 
supports the synchrotron nature of the nonther- 
mal emission in this remnant. 

Our fitted values for the break frequency are 
all in the range v c ~ (1.5 — 4) x 10 16 Hz. The en- 
ergy of an electron producing its peak synchrotron 
radiation at v in a magnetic field B is given by 

^ = 45 U7xlO"> Hz J (kT^gJ erg ' 

(1) 

Our values for E m then range from 14 to 45 erg (9 
to 28 TeV), at the low end of the range of upper 
limits to E m found for bright Galactic remnants 
by Reynolds & Keohane (1999). 

Since RCW 86 is a more highly evolved remnant 
than other objects in which synchrotron X-rays 
have been discovered, it is likely that the mecha- 
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Fig. 10. — ASCA spectrum of the Soft Region 
(SIS data - top, GIS data - bottom). A three- 
component model is also plotted (see Figure 9 for 
explanation) . 



Fig. 11. — ASCA spectrum of the Hard Region 
(SIS data - top, GIS data - bottom). A three- 
component model is also plotted (see Figure 9 for 
explanation) . 



nism causing the rolloff in the electron distribution 
above E m is synchrotron (and inverse-Compton) 
losses. Reynolds (1998) gives expressions for the 
maximum electron energy producible by shock ac- 
celeration limited by radiative losses, assuming 
that the electron mean free path is a factor rj times 
the gyroradius. In the strong-turbulence (Bohm) 
limit, r] ~ 1. For such conditions, with a shock 
compression ratio of 4, we have 

E m ~ 50 r 1 ' 2 (3^) u 8 erg (2) 

where B\ is the upstream magnetic-field strength 
and u% the shock velocity in units of 10 8 cm 
s _1 . Since optical inferences (Ghavamian 1999; 
Ghavamian et al. 2000) give nonradiative-shock 
velocities in the range of 400 to 900 km s _1 , we 
find predicted values for the loss-limited break en- 
ergy quite close to the break energies we infer from 
the fits, for Bi of a few microgauss and r\ ~ 1. 
A loss-limited spectrum will not have the exact 
shape of the srcut spectrum, but will roll off more 



slowly. However, given the limited spectral range 
of the fits, we cannot distinguish that difference, 
which in any case is less than the difference be- 
tween the srcut spectrum and a straight power- 
law. The coincidence between expected and ob- 
served cutoff energies, while not an ironclad con- 
clusion, adds circumstantial support to the idea 
that synchrotron X-rays are present in RCW 86. 
The electron-acceleration properties do not appear 
to vary much among the three regions we have 
studied; the primary difference in the nonthermal 
emission is just the overall strength of the emis- 
sion, perhaps as a result of differences in magnetic- 
field strength. 

The thermal X-ray emission consists of two dis- 
tinct components: a low-temperature component 
responsible for X-ray lines at low energies and a 
high-temperature plasma where Fe Ka line is pro- 
duced. The ratio of their emission measures sys- 
tematically varies from a high (~ 10) value in 
the Soft Region, through intermediate (4-5) val- 
ues in the Medium-Hard Region, and attains the 
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lowest (2) value in the Hard Region. The high- 
and low-temperature plasma should be in an ap- 
proximate pressure equilibrium, n e iT e i ~ n e hT e h, 
where n e i and n e h are electron densities in the 
low- and high-temperature components, and T e ; 
and T e h are their electron temperatures. For the 
Medium-Hard Region and the pshock+nei+sreut 
model, T el = 0.75 keV and T eh = 5.3 keV, and 
the assumption of a pressure balance gives a den- 
sity ratio n e i/n e h ~ 7. But their emission mea- 
sures differ by only a factor of 4, which suggests 
that the filling fraction fh of the high-temperature 
component is much higher that the filling fraction 
/; of the low-temperature component. If we as- 
sume that the hot gas fills most of the volume 
of the Medium-Hard Region, very approximately 
equal to 6 3 pc (where 6 pc is the linear size of the 
ASCA extraction region at the assumed distance 
of 3 kpc), then we arrive at n e h ~ 0.2 cm -3 (where 
fh was set to 1), based on the emission measure 
of the hot component listed in Table 2. Simi- 
larly, the low-temperature component has density 
n e i ~ 0.44/J^ 2 cm~ 3 , while we obtain n e i = 1.6 
cm~ 3 by assuming a pressure balance between 
the high- and low-temperature components. In 
order to make these estimates consistent, a low 
filling fraction of 0.075 is required for the low- 
temperature component. A low filling fraction 
fi is indeed suggested by a highly filamentary 
morphology of this region of the remnant seen in 
the archival HRI ROSAT images (Bocchino et al. 
2000) . The thickness of sharpest filaments seen in 
the HRI images is of the order of 10" or even less, 
which is much smaller than the 7' x 7' angular 
size of the Medium-Hard Region. 

The low-temperature X-ray component is most 
likely produced in the blast wave propagating into 
the ambient medium. The electron temperature 
T e i of 0.8 keV in the pshock+nei+sreut models 
implies a shock speed of 800 km s _1 if electrons 
are fully equilibrated with ions. In the absence of 
electron-ion equilibration, the shock speed might 
be even higher, above 1000 km s~ x as suggested by 
mean postshock temperatures in excess of 1 keV 
in the Sedov+nei+srcut models. The 800 km s _1 
shock velocity is at the upper velocity range of 
nonradiative, Balmcr dominated shocks seen in 
RCW 86 (Long & Blair 1990; Ghavamian 1999; 
Ghavamian et al. 2000). The preshock density 
n o = "ei/4-8 ~ 0.3 cm~ 3 is also in accord with 



estimates of the preshock density in fast nonradia- 
tive shocks in RCW 86. For example, Ghavamian 
(1999) finds preshock density of 1 cm -3 for a ~ 600 
km s _1 shock located at the southern boundary of 
the Hard and Medium-Hard Regions. The shock 
age, t <~ Ti/n e i — 3, 600 yr, is in reasonable agree- 
ment with the remnant age estimated at ~ 10 4 
yr by Rosado et al. (1996) from standard argu- 
ments based on the Sedov dynamics, and at 4,300 
yr by Petruk (1999) whose estimates were based 
on 2-D hydrodynamical simulations. This agree- 
ment further supports the identification of the low- 
temperature thermal component with the material 
shocked by the blast wave. Bocchino et al. (2000) 
also attributed the low-temperature component to 
the shocked ambient gas based on analysis of Bep- 
poSAX and ROSAT data. 

The high-tcmperature thermal component is 
produced by hot, tenuous gas in the remnant's 
interior. The best evidence for this is provided 
by the BeppoSAX hardness map (Bocchino et al. 
2000), which shows that the hard X-ray emission 
is located interior to the soft X-ray emission. The 
high (5.3 keV) temperature of the hot component 
suggests that this material was thermalized in fast, 
several thousand km s" 1 shocks at earlier stages 
of evolution of RCW 86. As determined by fit- 
ting GIS data in § 2, its ionization age is less 
than 6.7 x 10 9 cm -3 s at the 90% confidence level, 
which in combination with its electron density of 
n e h = 0.2 cm~ 3 leads to an upper limit of 1000 
yr for its age. However, if Fe resided in dust 
grains instead of gas, this upper limit would be 
several times higher, and entirely consistent the 
blast wave age derived previously. The Fe abun- 
dance in the hot gas is uncertain, but may be en- 
hanced with respect to the solar Fe abundance. 
If real, this would imply the presence of heavy- 
element SN ejecta in the hot thermal component, 
as already suggested by Bocchino et al. (2000). 
This component might have originated in the in- 
teraction of the blast wave with a low-density cav- 
ity created by the SN progenitor prior to its ex- 
plosion (Vink et al. 1997). It cannot be the nor- 
mal hot interior of a Sedov blast wave, at least 
for the Sedov profile corresponding to the current 
blast- wave conditions, or our previous Sedov fit- 
ting would have accounted for it. 

Because of the simultaneous presence of two (or 
more) thermal components and of the nonthermal 
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continuum, the separation of the X-ray spectrum 
of RCW 86 into individual components is difficult 
and uncertain. The separation of thermal emis- 
sion into a low-temperature single shock compo- 
nent and a single ionization timescale component 
with high (and constant) temperature is certainly 
a reasonable first approximation, but one might 
expect that such a simple model would not be 
able to provide a perfect match to the observed 
X-ray spectrum. This is indeed the case because 
our best fits are not statistically acceptable at low 
photon energies, where ASCA provided a high- 
quality X-ray spectrum with the excellent S /N ra- 
tio. We speculate that these problems stem from 
the simultaneous presence of shocks with differ- 
ent velocities in our extraction regions. Optical 
observations indicate presence of shocks in a wide 
range of shock velocities from ~ 100 km s^ 1 to 900 
km s _1 in various regions of RCW 86. Both slow 
(~ 100 km s- 1 ) and fast (~ 800 km s" 1 ) shocks 
are within the extraction regions, as evidenced by 
the presence of radiative optical filaments and by 
the ASCA spectra, respectively. While ~ 100 
km s _1 shocks cannot produce X-ray emission, 
intermediate- velocity (~ 600 km s _1 ) shocks ca- 
pable of producing X-ray emission are present at 
the southern boundary of the Medium-Hard and 
Hard Regions (Long & Blair 1990; Ghavamian 
1999; Ghavamian et al. 2000). X-ray emission 
from such intermediate- velocity shocks might have 
been already detected in the ROSAT PSPC spec- 
tra of RCW 86 by Bocchino et al. (2000), who 
reported the presence of very low (0.29 keV) tem- 
perature gas in the SW region of RCW 86. X-ray 
emission from these intermediate-velocity shocks 
is not accounted for in our simple decomposition 
of the X-ray spectra. Our derived model parame- 
ters must then be considered as preliminary rather 
than as robust estimates. In particular, low post- 
shock electron temperatures in Sedov models (Ta- 
ble 2) might be an artifact of the model used 
instead of a measurement of the postshock elec- 
tron temperature. In view of the limited spatial 
and spectral resolution of the ASCA satellite and 
its detectors, it would be premature at this time 
to construct a more sophisticated thermal models 
with multiple shock components. 



4. DISCUSSION 

Detailed modeling of X-ray spectrum of RCW 
86 revealed the presence of the nonthermal syn- 
chrotron emission, in addition to high- and low- 
temperature thermal components. This nonther- 
mal emission is strongly correlated with the radio 
emission, and the derived break frequency v c of 
2 — 4 x 10 16 Hz is consistent with the maximum 
electron energy producible by shock acceleration 
limited by radiative losses. The presence of non- 
thermal synchrotron emission explains why previ- 
ous attempts to model X-ray spectrum of RCW 86 
in terms of pure thermal emission led to puzzling 
results. 

The nonthermal synchrotron emission is most 
important at intermediate (1-3 keV) photon ener- 
gies, while at low (high) energies the ASCA X-ray 
spectrum is dominated by a low- (high-) tempera- 
ture thermal component. This is a more complex 
behavior than seen in remnants such as SN1006 
where thermal emission dominates at low energies, 
while nonthermal emission is most pronounced at 
high energies. To model the nonthermal compo- 
nent, a physically reasonable description is neces- 
sary, as much for the implied inferences about the 
thermal emission as for the synchrotron emission 
itself. We do not expect power-laws in X-ray syn- 
chrotron emission as we are always observing the 
high-energy rolloff; using power-laws for modeling, 
while perhaps allowed by the data in some cases, 
can lead to misleading inferences about spectral- 
line equivalent widths. For the srcut model with 
a break frequency of 3.7 x 10 16 Hz, a power-law 
tangent to that model at 2.1 keV (in the region 
where synchrotron emission dominates in RCW 
86) overpredicts the continuum at 6.4 keV by 50% 
compared to srcut - so the equivalent width of the 
Fe Ka line would be considerably underestimated 
by a power-law fit. 

The nonthermal spectrum traverses the IR and 
optical on its way from radio to X-ray wave- 
lengths. Since all the fitted break frequencies are 
well above the frequency of the center of the V 
band (5.5 x 10 14 Hz), we can assume a straight 
extrapolation from radio frequencies with a slope 
of 0.6. If the 4.5 Jy at 1 GHz we attribute to 
the 7' x 7' Medium-Hard region were spread uni- 
formly, they would imply a surface brightness at 
V of about 29 m arcsec -2 . However, the morphol- 
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ogy should be identical to the radio; in particular, 
the radio emission is concentrated in quite thin fil- 
aments, so that the peak brightness is larger than 
average by at least a factor of 10, raising the pre- 
diction to 26 m .5 arcsec -2 , perhaps attainable. At 
2 microns, surface brightness would be almost a 
magnitude brighter. 

The derived blast wave velocity of 800 km s^ 1 , 
the estimated preshock density of ~ 0.3 cm~ 3 , and 
the shock age of ~ 4 x 10 3 yr are consistent with 
previous estimates, although a low (0.075) filling 
fraction for the gas shocked by the blast wave 
seems to indicate that emission is dominated by 
denser than average material. The enhanced X- 
ray emission and the presence of radiative shocks 
in the bright SW section of the remnant have been 
attributed to the interaction of the blast wave with 
a dense (~ 10 cm~ 3 ) interstellar cloud. Our find- 
ings support this interpretation. The presence of a 
very hot (~ 5 kcV) gas with a large filling fraction, 
with the emission measure far in excess of that ex- 
pected in a standard Scdov model, indicates that 
the SN progenitor most likely exploded within a 
low-density cavity. This suggests that RCW 86 is 
a remnant of a core-collapse SN, whose progenitor 
was perhaps a member of an OB stellar associa- 
tion seen at a distance of 2.5 kpc in this region of 
the sky (Westerlund 1969). 

The presence of nonthermal synchrotron emis- 
sion in RCW 86 is certainly one of the most inter- 
esting aspects of this well-studied remnant. Ener- 
getic electrons producing this emission have been 
presumably accelerated in the blast wave, so it 
would be of great interest to study in greater de- 
tail the relationship between the thermal and non- 
thermal emission, and between optically-emitting 
shocks and the nonthermal emission. The correla- 
tion of both thermal and nonthermal X-ray emis- 
sion with the strong far-IR emission from collision- 
ally heated dust grains seen in RCW 86 (Dwek et 
al. 1987; Arendt 1989; Greidanus & Strom 1990; 
Saken et al. 1992) is also important in view of the 
crucial role of dust grains in acceleration of cos- 
mic rays (Meyer, Drury, & Ellison 1997; Ellison, 
Drury, & Meyer 1997). Such studies would allow 
for better understanding of the shock acceleration 
process. Unfortunately, the quality of the present 
X-ray observations does not allow for a unique 
separation between the thermal and nonthermal 
emission. In order to achieve this separation, fu- 



ture observations of RCW 86 with a spatial and 
spectral resolution higher than that of ASCA are 
necessary, for example with the XMM-Newton or 
Chandra satellites. ~XMM.-Newton should provide 
high S/N CCD spectra with a moderate spatial 
resolution, while Chandra might be able to resolve 
even the most compact X-ray filaments seen in the 
ROSAT HRI images. The brightest X-ray filament 
seen in the ROSAT HRI image of the SW region of 
RCW 86 (Bocchino et al. 2000) is of particular in- 
terest, because it coincides with the brightest non- 
radiative shock in RCW 86 whose shock velocity 
and the preshock density was already determined 
from the optical spectroscopy (Ghavamian 1999). 
It is also critical to obtain higher spectral reso- 
lution data, at the resolution provided by the X- 
ray microcalorimeter onboard of the failed Astro-E 
satellite. Such observations would allow for deter- 
mination of gas temperature from line diagnostics, 
and hence for unambiguous separation of thermal 
and nonthermal components. Once these compo- 
nents are separated, one can combine information 
from both components to learn about acceleration 
of energetic (TeV) electrons in SNRs. Maybe we 
can then find out why the nonthermal contribu- 
tion is so important in RCW 86, and relatively 
weak in some SNRs. 
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